Correlation effects in electronic structure of PuCoGa 5 
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We report on results of the first realistic electronic structure calculations of the Pu-based PuCoGas 
superconductor based on the dynamical mean field theory. We find that dynamical correlations 
due to the local Coulomb interaction between Pu /-electrons lead to substantial modification of 
the electronic structure with a narrow peak being formed in vicinity of Ef, in agreement with 
the experimental photoemission spectra, and in contrast with the recent calculations within the 
LDA+U method y|, where only static electronic correlations have been included. Both Pu and Co 
contribute in equal footing to the narrow peak on the density of states at the Fermi level, the Co 
partial density of states being prominently affected by electronic correlations on the Pu sites. The 
k-resolved spectral density is calculated and the theoretical spectral function resolved extended Van 
Hove singularity near Ef- This singularity may lead to enchancement of the magnetic susceptebility 
and favour d— wave superconductivity. 
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The recent discovery of superconductivity in 
PuCoGas and PuRhGa 5 Q has attracted a lot of 
attention to electronic structure of these Pu-based 
compounds. They have the same HoCoGas-type 
tetragonal structure as the Ce-based superconduc- 
tors CeMGas (M is transition metal) however 
PuCoGas has an order of magnitude higher transi- 
tion temperature T c =18.5 K, the highest T c among 
known /-electron superconducting materials. Both 
CeAfGas and PuAfGas are unconventional super- 
conductors with the nodal d— wave symmetry of 
superconducting gap A and antiferromagnetic fluc- 
tuations being a likely mechanism for their uncon- 
ventional superconductivity. It has been suggested 
|(| that the increase of T c in the PuMGas com- 
pounds in comparison with CeMGas might be due 
to more itinerant character of the 5/ band in Pu. 

The experimental data show no evidence of a 
long-range magnetic structure down to temperature 
of 1 K, whereas the magnetic susceptibility obeys 
the modified Curie- Weiss law with an effective lo- 
cal moment 0.68 The specific heat coefficient 
7=77 mJ mol _1 K -2 is slightly larger than the 
value observed in <5-Pu (50 mJ mol -1 K -2 0) in- 
dicating possible enhancement by strong dynamical 
spin fluctuations. The experimental photoemission 
spectra (PES) of PuCoGas Q exhibits a broad fea- 
ture centered 1.2 eV below the Fermi level, which 
is made up of Co d- and Pu /-bands, as well as a 
narrow peak at the Fermi level. The latter feature is 
similar to narrow resonances observed in the PES of 
<5-Pu as well as in the spectra of some Pu-based 
compounds fToL fTT| - 

The electronic structure of PuCoGas has been in- 
vestigated in several papers on the basis of the local 
spin density approximation (LSDA) [Tl 111 HI HH . 
In the non-s pin- polarized band structure calcula- 
tions of Refs. [13, E f° r nonmagnetic state the main 



/-state spectral weight is located in the vicinity of 
the Fermi level in contradiction with experiment. 
The spin polarization pushes occupied part of the 
/-band below the Fermi level [LH [l4j , with the den- 
sity of states (DOS) for the antiferromagnetic phase 
1 1 ll ] being in good agreement with the experimental 
PES 0. However, the Pu total magnetic moment 
obtained within the LSDA is too large (around 2 jug) 
in comparison with experimentally measured fi e f f as 
well as the ordered magnetic structure in PuCoGas 
predicted by the LSDA calculations has not been ob- 
served in experiment. Joyce et al. [8j obtained theo- 
retical PES in a very good agreement with the exper- 
imental one using the mixed level model (MLM) . In 
the MLM 5/ electrons of Pu is divided in four local- 
ized and one itinerant, however the energy position 
of the localized states is not calculated ah initio but 
rather treated as a parameter which value is chosen 
to match experiment. 

The strong enhancement of the specific heat coef- 
ficient 7, as well as the characteristic narrow reso- 
nance at the Fermi level, hint on possible importance 
of dynamical electronic correlations in PuCoGas. In 
<5-Pu 0| and Pu monochalcogenides 0] a good 
agreement between the theoretical and measured 
PES for non-magnetic phase (which is the ground 
state for these materials in accordance with the ex- 
periment) has been obtained only when dynamical 
correlations between /-electrons were taken into ac- 
count by means of the dynamical mean-field theory 
(DMFT) [ill in a framework of the LDA+DMFT 
[HIH method. For PuCoGa 5 Shick et al. have 
taken into account only static correlations by means 
of the LDA+U method. They obtained the antifer- 
romagnetic ground state with the reduced (in com- 
parison with the LSDA) Pu effective magnetic mo- 
ment around 1 (13- At the same time, the density 
of states reported in Ref. Q demonstrates some dis- 
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crepancy with the experimental PES. While the po- 
sition of the broad manifold reproduced quite well, 
the narrow feature at Ep is absent. It seems reason- 
able to suggest that by analogy with (5-Pu and other 
Pu-based compounds an adequate description of the 
electronic structure of PuCoGas requires dynamical 
correlations to be taken into account. It is impor- 
tant to notice that in the Pu-based compounds the 
dynamical correlations affect most profoundly the 
states in vicinity of the Fermi level |17| , which are in- 
volved in superconductivity. It is the purpose of this 
Letter to present electronic structure of PuCoGas 
with the dynamical electronic correlations explicitly 
included within the LDA+DMFT approach. 

We start with the LSDA+U calculations of 
PuCoGas. Following Ref. Q to allow for the anti- 
ferromagnetic ground state we double the unit cell 
along the a and b crystallographic axes and point 
magnetic moments of the nearest- neighbor Pu atoms 
in the Pu planes in opposite directions, the mag- 
netic moments of Pu atoms being aligned along the 
crystallographic c axis. The calculations have been 
carried out by the full-potential linear MT-orbitals 
method (FPLMTO) [2l| at experimental lattice pa- 
rameter of PuCoGas (a=7.845 A, c/a=1.603 for the 
unit cell) with the spin-orbit coupling included in 
a self-consistent second-order variational procedure. 
We include the full matrix of on-site Coulomb in- 
teraction between Pu 5/-electrons parametrized by 
the average Coulomb interactions U =3 eV and the 
Hunds-rule exchange J ==0.55 cV. These values of 
U and J are in the range of commonly accepted val- 
ues for Pu. For the double counting (DC) term we 
have employed so-called "around mean-field" formu- 
lation of the LDA+U method [H This ver- 
sion of the LDA+U method was recently successfully 
used to explain non-magnetic ground state of the Pu 
5-phase [24| . We used previously the same values of 
U and J as well as the same choice of the DC in the 
DMFT simulations of Pu monochalcogenides [lTj . 

The one-particle Hamiltonian H t (k) obtained by 
the converged LDA+U calculations was then orthog- 
onalized by the Lowdin transformation |25| and used 
to compute the local Green's function by means of 
the Brillouin zone (BZ) integration 

G(iu>) = + - fft(k) - £ dc Mr\ (1) 

k 

where u> = (2n + l)irT are the fermionic Matsubara 
frequencies for a given temperature T, fj, is the chem- 
ical potential and T, dc (iu>) is the local self-energy 
with a "double counting" term subtracted. Within 
the DMFT scheme the local self-energy has been ob- 
tained by the solution of the many-body problem for 
a single quantum impurity coupled to an effective 
electronic "bath" through the Weiss field function 

m 

g-\iuj) = G~ x {iu) + H dc {iLu). (2) 
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FIG. 1: The total DOS (solid line) as well as the par- 
tial DOS for Pu (dashed line) and Co (dot-dashed line) 
atoms obtained within the DMFT. The thick solid and 
dot-dashed lines are the experimental >8j and theoretical 
PES in arbitrary units, respectively. 

As a quantum impurity solver we have employed 
the recently developed spin-orbit version [Tg of 
the T-matrix and fluctuating exchange technique 
(SPTF) [HI23. In the spin-orbit SPTF the local 
self-energy is written as a sum of three contributions: 

S = E (TH) +E (TF) +5] (PH) 5 (3) 

where £< Tff ) and £< TF ) are the T-matrix "Hartree" 
and "Fock" contributions, respectively, Y,( PH ^ is the 
particle-hole contribution. E( Ti? ) and £( TF ) are ob- 
tained by substitution in the corresponding Bartree 
and Fock diagrams the bare Coulomb interaction 
with the frequency dependent T-matrix, the latter 
is given by summation of the ladder diagrams in the 
particle-particle channel. ^ PH ^ is obtained by the 
RPA-type summation in the particle-hole channel 
with the bare vertex being substituted by static limit 
of the T-matrix 0, [^3 • While being of perturba- 
tive type the spin-orbit SPTF solver has been shown 
to provide an adequate description of correlations in 
3d metals 26, 27], half-metallic ferromagets [2^], and 
some actinide systems [l7| . 

For the double counting term in DMFT simula- 
tions we have used the static limit of the self en- 
ergy For the spin-orbit SPTF quantum impu- 
rity solver we employed 10 3 Matsubara frequencies 
and temperature 500 K. It is important to notice 
that with our choice of the double counting term 
the Bartree-Fock contribution is already included 
into the LDA+U one-particle Bamiltonian, therefore 
making magnetization almost temperature indepen- 
dent. We have carried out DMTF iterations until 
convergence in both the chemical potential jj. and 
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FIG. 2: The spectral function of PuCoGas. The sequence of colors in order of increasing spectral density is blue, 
green, yellow, and red. In the inset the BZ of PuCoGas is shown together with our designations for the high-symmetry 
points. 



the local self-energy was achieved. Then the Pade 
approximant is used p9j for analytical continuation 
of the local self-energy to the real axis in order to ob- 
tain the DOS n(E) = -1/ttTtQ[G(E)} and spectral 
function A{E, k) = -l/>rrTr3[G(E, k)]. 

We started from the antiferromagnetic state of 
PuCoGas but the self-consistent calculations con- 
verged to practically zero magnetic moments both 
on the Pu and Co sites. The long-range magnetic 
order is suppressed due to combined influence of a 
strong rclativistic spin-orbit coupling and the local 
Coulomb interaction effects. The j^ji states of Pu in 
PuCoGas are almost completely filled (the occupa- 
tion number for /-band n^=5.5) whereas the empty 
/ 7 / 2 states are located around 3 eV above the Fermi 
level. It is this splitting that suppress the magnetism 
of Pu in comparison with the LSDA electronic struc- 
ture. The Co 3c? band is located between -2 eV and 
the Fermi energy. 

The DMFT results for PuCoGa 5 DOS is shown 
on the Fig. ^ together with the experimental PES 
spectrum || and the theoretical PES obtained from 
the DMFT DOS. In order to compute the theoretical 
PES spectrum the calculated DOS data have been 
multiplied by the Fermi function corresponding to 
the experimental temperature of 77 K and then 
convoluted with a Gaussian to simulate for the in- 
strumental resolution and lifetime broadening of the 



core hole |g, y(j ■ One may see a pronounced quasi- 
particle peak at the Fermi level in the DMFT DOS. 
The theoretical PES agrees quite well with the ex- 
perimental one, the shape of the feature at the Fermi 
level as well as relative spectral weight of the broad 
manifold and the resonance at Ep being reproduced 
very accurately. However, the manifold is somewhat 
wider in the theoretical PES compare with the ex- 
periment, with a peak at -0.6 eV and a shoulder at 
-1.8 eV. 

Interestingly, while in our DMFT calculations the 
correlation effects have been taken into account only 
on the Pu sites, both the Pu /- and Co d-bands 
are substantially modified in comparison with the 
LDA+U result and they both contribute in equal 
footing to the narrowpeak at Ep. The LSDA [15j 
as well as LSDA+U [l| calculations for the antifer- 
romagnetic phase predict that the Pu /-band gives 
by far the dominant contribution to DOS at Ep. 
We have performed also the LDA+DMFT calcula- 
tions with the local Coulomb interaction between 
Co d-electrons being included as well (with U =2 eV 
and J=0.9 eV), however that leads to rather minor 
changes in the electronic structure. 

The spectral function of PuCoGas calculated 
within the DMFT is displayed in Fig. The max- 
imums in spectral function around Al, A and Z 
points correspond to both the Pu /-band and Co 
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d-band, contributing mainly to the manifold at en- 
ergies between -2 and -1 eV. The increase of the spec- 
tral weight near the the Fermi level due to van Hove 
singularities (saddle points) close to R and T points 
is mainly due to the Co d-band, and in between Z 
and R points is due to contributions from both Co 
and Pu. The flat Co d band located exactly at the 
Fermi level is clearly seen in the spectral function 
along the R — X direction. The extended van Hove 
singularity near R point together with the van Hove 
singularity at T point result in strong q-dependence 
of magnetic susceptibility (cf., e.g., the "parquet" re- 
sults for the two-dimensional Hubbard model [3l|p. 
This enhanced susceptibility can lead to d— wave su- 
perconductivity [33|. Relevance of the van Hove 
singularities for the superconductivit y ha s been dis- 
cussed for high-temperature cuprates |3j|33, as well 
as for Sr2Ru04 |34j . However, in comparison with 
these cases, the electron spectrum of PuCoGas is a 
bit less anisotropic and probably can not be con- 
sidered as a purely two-dimensional one even in a 
zeroth-order approximation (see the results for R— X 
direction in Fig. El> ■ 

In conclusion, we have shown that the dynam- 
ical correlation effects modify electronic structure 



of the PuCoGas superconductor in essential way. 
We have taken the dynamical correlations into ac- 
count by means of the dynamical mean-field the- 
ory using the spin-orbit T-matrix and fluctuating 
exchange (SPTF) approximation for the quantum 
impurity solver, in agreement with the experimen- 
tal data, the self-consistent solution turned out to 
be non-magnetic. The DOS obtained by the DMFT 
simulations is in a good agreement with the exper- 
imental PES, with the narrow peak at the Fermi 
level being very well reproduced. While only Pu 
sites have been treated as correlated, the dynamical 
correlations modify the Co bands as well. On the 
other hand, taking into account local correlations on 
Co sites does not effect noticeably on the electronic 
structure. Both the Co and Pu bands contribute al- 
most on equal footing to the quasiparticle resonance 
at Ep, in contrast with the LSDA picture, where 
only Pu bands contribute at DOS on the Fermi level. 
Theoretical spectral function demonstrates the ex- 
tended Van Hove singularity near R point which is 
due to the Co d band. New ARPES experiments are 
highly desirable to investigate the energy spectrum 
of PuCoGas in the vicinity of the Fermi level and 
check the theoretical predictions. 
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